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1. Introduction 

Since its formulation IIBI. the Bohr-Mottelson model remained the indispensable tool for 
interpretation of the collective phenomena in even-even nuclei, with numerous solutions 
being constantly proposed by considering different shapes and approximations for the 
potential energy V{/3,'y) [3H5]. These simple solutions generated in turn other more 
involved approaches by considering octupole degrees of freedom [BHS] and deformation- 
dependent mass terms [9l[Tn], or more simple ones by imposing 7 -rigidity [UllIS]. The 
formalism of the Bohr-Mottelson model recently inspired similar approaches for the 
dehnition of collective Hamiltonians related to such phenomena like the chiral modes [13] 
or the wobbling excitations m. The reformulation of the Bohr Hamiltonian in six 
dimensions [I5] is another recent noteworthy addition to the geometrical description of 
collective motion. 



Conjunction of'j-rigid and '^-stable collective motion... 


2 


The parameter independent solutions of the Bohr Hamiltonian are especially 
useful in the sense that they serve as reference points for the energy spectra and 

electromagnetic properties of nuclei. It was found that this is also the case for the 

critical point symmetries E{b) [16] and X(5) [I^ describing the shape phase transitions 
from vibrational to 7 -unstable and respectively to the axially deformed nuclei. The E{b) 
solution is exact, while X(5) employs two approximations, one related to the separation 

of variables and the other based on the small angles for the 7 shape variable. The 

domains of validity for these approximations were found to be incongruous when the 
Hamiltonian is numerically diagonalized [18] in the framework of the recently developed 
Algebraic Collective Model [I9H2I]. Also, due to its structure, -A(5) is parameter free 
only for the ground and /3 bands, while the 7 band is determined by the stiffness of the 
7 oscillations. All these shortcomings were commented on all accounts in [ 22 ], where 
an exactly separable version of the X(5) model denoted ES-X(5), was proposed. As in 
this case the /3 and 7 degrees of freedom are treated on the same footing, the model 
will no longer be parameter independent, the whole energy spectrum depending on the 
7 stiffness parameter. The 7 -rigid version of this model, called X(3) [12], however will 
be parameter independent. Although the 7 -rigidity hypothesis is somehow a crude one, 
it provides simple approaches with successful reproduction of the relevant experimental 
data [II11I21123H27]. Moreover, an unexpected similarity between the (3 excited bands 
of X (5) and X (3) solutions was found, which inevitably address the question about the 
importance of rigidity in explaining the critical collective phenomena. 

In the present work we propose a simple exactly separable model for the competition 
between the 7 -rigid and 7 -stable collective motion in the phase transition between 
spherical and deformed shapes. The coupling of the two types of (3 vibration is achieved 
by introducing a control parameter measuring the degree of the system’s 7 -rigidity in 
an Ising type Hamiltonian [1]. The separation of variables is achieved by considering 
a potential of the form u{/3) + u{'y)//3^ [3112 ^1281130] adapted to the current problem. 
Matching the two competing excitations, the 7 potential is chosen to be a harmonic 
oscillator centered in 7 = 0 , which is consistent with the prolate 7 -rigid part of the 
problem. As the purpose of this paper is to investigate the 7 rigidity near a critical 
point which is commonly identihed with a very flat potential, u{/3) is taken to be an 
inhnite square well. Although the separation of the potential leads to independent 
analytical descriptions of the (3 and 7 degrees of freedom, it contravenes the Liquid 
Drop Model prescription for the potential energy However, the chosen potential 
shape can be simulated by considering higher order anharmonic terms in the generalized 
collective geometrical model [32] . 

The energy spectrum and E2 transition probabilities resulting from the proposed 
model depend up to an overall scaling factor on two parameters, namely the rigidity and 
the stiffness of the 7 vibrations. Their separate influence on the model’s characteristics is 
investigated through numerical applications. The experimental realization of the model 
is found in few rare earth nuclei around N = 96. 
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The kinetic energy operator associated to a prolate 7-rigid nucleus is 
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while in the 7-stable case it is expressed as 
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where Q is the angular momentum operator from the intrinsic frame of reference with 
the corresponding components Qk{k = 1 , 2 , 3 ), while B is the mass parameter. The 
interplay between 7-stable and 7-rigid collective motion is achieved by considering the 
Hamiltonian: 

H = xTr + {l-x)Ts + V{/d,P, ( 2 . 3 ) 

where 0 ^ y < 1 is a parameter which indicates the degree of the system’s rigidity 
against 7 vibrations. Note that the prolate 7-rigid limit y = 1 is avoided in order to 
preserve the geometry of the curvilinear space. 

In order to achieve an exact separation of the fd variable from the 7-angular ones, 
the reduced potential is considered of the form similar to that used in [ 3 ll 2 ^[ 28 l - [ 3 n] but 
adapted for the present problem: 


<’(/?. 7 ) = ^r(/3,7) = ti(/3) + (i-x)'pr 


( 2 . 4 ) 


Factorizing the total wave function as T(/3,7,r2) = ^(/3)(p(7, hi), the associated 
Schrodinger equation is separated in two parts: 
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where e = ^E. As the present model is based on the hypothesis of small oscillations 
around 7 = 0 which cease at y = 1, the rotational term of the last equation is 
approximated as in HZI: 
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Contrary to the 7-unstable case [ 33 ] , the 7 and angular variables can also be separated 
by considering the product state (^(7,12) = Here are the 
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Wigner functions associated with the total angular momentum L and its projections 
on the body-hxed and laboratory-hxed 2 ; axis, M and K, respectively. Averaging the 
approximated equation (12.61) on this product state one obtains the following equation 
for the 7 shape variable: 

h(7) = evil), (2. 

with 


Id d 
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sm 37 ay ay 4 sm y 
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(2.9) 


For small oscillations around y = 0, the ’’frozen” value of the y-rigid part of the 
Hamiltonian, it is natural to adopt a harmonic oscillator form for the y potential, 

“(7) = 


where the parameter a dehning the string constant of the oscillator is usually referred 
to as the stiffness of the y vibrations. Then fl2.8p can be treated as in [HI[22], i.e. 
by applying a harmonic approximation for the trigonometric functions around y = 0 . 
This leads to a differential equation for y which resembles the radial equation for a two 
dimensional harmonic oscillator: 
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The solutions are readily obtained in terms of the Laguerre polynomials nil: 

( 2 . 12 ) 

where Nn,\K\ is a normalization constant and n = {n.y — \K/2\)/2. The corresponding 
eigenvalues are 


e.y = + 1), n.y = 0,1 , 2 ,..., (2.13) 

with K = 0, ±2n.y for n.y even and K = -ii2n^ for odd, respectively. At this point it 
would be opportune to comment the implications of choosing such a y potential. First 
of all it loses the symmetry properties required by the Liquid Drop Model prescription 
for the potential energy [1]. This results from the fact that the harmonic oscillator 
potential (I2.10p is obtained from the second order approximation of the symmetry 
obeying potential a^(l — cos3y). Secondly, the corresponding Hamiltonian is no longer 
hermitical with respect to the scalar product dehned with |sin3y|(iy measure |34j . 
However, the adopted approximation was demonstrated in [TH| to be reliable for high 
stiffness, condition which is achieved in the present model as will be shown in the 
numerical application. 

In what concerns the /3 degree of freedom, one will consider here an anharmonic 
behaviour reflected into a square well shape of the potential: 
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where fdw indicates the position of the infinite wall. In this way, by varying the rigidity 
parameter from y = 0 to y —)■ 1, the system will change from ES-X(5) model | 22 ] to a 
qnasi-y-rigid solution whose y-rigid limit is the X (3) model [12]. Now, equation (12.hh can 
be brought to a Bessel differential equation by the change of variable ^{(3) = 


Id , ,, . 


a /32 




fm = 0 , 


where 
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L(L + 1) - (1 - x)K'‘ ^ ^ 
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The boundary condition f{/3w) = 0 gives the (3 energy spectrum in terms of the s-th 
zero Xs,u of the Bessel function Ju{xs,uf3/f3w)'- 
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The order of the Bessel function’s zero is related to the (3 vibration quantum number by 
Up = s — 1. The excitation energy of the whole system in respect to the ground state is 
then defined as 




EK,np,n-,{L) — — [eL,K,np+i,n.y{f3w) “ eo,o,i,o(/3m)] , 


(2.18) 


where the set of quantum numbers {K,np,n-y} uniquely identify a rotational band, i.e. 
{000} for the ground band, {010} and {201} for the hrst (3 and 7 excited bands and so 
on. Correspondingly the (3 variable normalized wave function is then given as: 

^L,K,np,n.yif3) = Nnf^,uf3^ ^ Ju{Xnp+l,ul3 / f3w), (2.19) 

where is the normalization constant obtained from the condition 
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With this fuuctiou oue cau calculate the trausitiou rates, employiug the geueral 
expressiou for the quadrupole trausitiou operator. 


Tf 2) = tfj 


1 


cosy + ^ {Dl^{n) + Dl_^{n)) siuy 


( 2 . 22 ) 


where t is a scaling factor. The hnal result for the E2 trausitiou probability is giveu in 
a factorized form 
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where C is the Clebsch-Gordan coefficient dictating the angular momentum selection 
rules, while B and G are integrals over the shape variables [3 and 7 with integration 
measures and respectively IsinSyl dy. Note that in case of the /S.K = 0 transitions 
described by the hrst term of E2 operator fl2.22p . the 7 integral reduces to the 
orthogonality condition for 7 n^,|x|( 7 ) wave functions. 




Figure 1. i? 4/2 (a), /3 (b) and 7 (c) band heads normalized to the energy of the first 
excited state are given as function of the rigidity parameter y for few values of the 
stiffness parameter a, namely 5, 10, 20 and 40. 

3. Numerical application and discussion 

The proposed model have two adjustable parameters not counting the scale, namely the 
stiffness a of the 7 vibrations and the rigidity y. Indeed, considering the energy spectrum 
and the E2 transition probabilities normalized to the energy of the hrst excited state 
and respectively to its B{E2) rate to the ground state, the dependence on the mass B 
and the position of the inhnite wall f3w vanishes. The evolution as function of y and 
a of theoretically evaluated spectral observables such as i? 4/2 = E{4C)/E{2'^) ratio and 
the (3 and 7 band heads normalized to the energy of the hrst excited state is visualised 
schematically in hgure [H From this hgure one can see that the hrst two sets of curves 
corresponding to i? 4/2 and E{Qi'^)/E{2'^) collapse at y —)■ 1 to the values of the X(3) 
energy spectrum, namely 2.44 and respectively 2.87. Although the ^(3) solution by 
construction excludes the existence of the 7 band, our approach at y ^ 1 provides the 
value 1 for the 7 band head. At y = 0, the corresponding values of ES-X(5) model 
are recovered for all three ratios. However, up to the y —)■ 1 limiting case the ratio 
i? 4/2 as well as the two excited band heads have quite diherent behaviours. Indeed, the 
functions are not linear but becoming more straighter for /5, and even more for the 7 
band head. Another diherence between the behaviour of the R 4./2 ratio and that of the 
excited band heads is the distinct effect of the stiffness parameter a. The increase in 
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stiffness for i? 4/2 is reflected in a longer and more horizontal plateau defined up to the 
turning point above which the line goes abruptly to the X(3) value. As a matter of 
fact the height of the plateau seems to reach a certain saturation at higher values of the 
stiffness, the saturation maximal height being Ri /2 = 3.33 which is the SU{3) rotational 
limit. In contradistinction, the effect of the stiffness increase on the excited band heads 
produce an almost proportional increase for the corresponding function. 

In what concerns the model’s wave functions, these are identical with those of ES- 
X(5) and X(3) when x = 0 and y ^ 1, respectively. But in the y ^ 1 case the 
norm of the /3 state is not the same as in the ^(3) model due to the different metrics 
of the deformation space, i.e. the integration measure in the present case is fd^dlS in 
comparison to fd'^dld of the X(3) formalism. As a consequence of this fact, a reshaping of 
the {3 probability distribution takes place. Indeed, as can be seen from the upper panels 
of hgure O while the probability peak in the X (3) case is positioned in the center of the 
square well, in the present model it is shifted to higher deformation and also becomes 
more pronounced. From the lower panels of hgure [2] one observes that this behaviour 
is perpetuated when going to smaller values of parameter y, reaching at some point a 
stagnation in the position and hight of the peak which continues up to the ES-X(5) 
limit. Increasing the stiffness, the stagnation interval diminishes and the peak at y = 0 
is shifted to a higher f3 deformation. This is the result of the centrifugal contribution 
from the 7 dependent term of the potential (I2.4p which compresses the wave function 
to higher (3 values as the stiffness of the 7 vibration increases. Such a behaviour of the 
(3 wave function is very much similar to the /3-rigid case na. In conclusion, the price to 
be paid when transforming 7 from a simple parameter to a variable is the conhnement 
of the ground state (3 distribution probability at higher expectation values of the [3 
deformation. 

An extensive search for nuclei which fall in the category described by the present 
model was made by means of model hts on rare earth and actinide isotopes with 
-R4/2 ^ 2.44 and with all three bands experimentally accounted. The quality of the 
hts was judged by the quantity 


a = 



77 ^( 2 +) 




(3.24) 


where E™ is the theoretical energy dehned by (I2.18p with k indexing angular momentum 
states from the ground, (3 and 7 bands. The experimental realization of our model was 
found to occur in ^®°Gd, ^®^Dy and ^®®Er isotopes. As can be deduced from table [U 
where the results of the hts are given in comparison to the available experimental data, 
the best agreement is obtained for ^®°Gd and ^®^Dy nuclei. For ^®^Dy one obtained 
the highest degree of mixing between 7 -rigid and 7 -stable rot at ion-vibration, judging 
by the rigidity value y = 0.269 which is associated with a comparable weighting for 
the two competing excitations. The hts for the other two nuclei provide much higher 
values for y, corresponding to prevalent 7 -rigid behaviour. Especially high 7 -rigidity is 
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Figure 2. The ground state /3 probability density in respect to the (i/3 integration 
measure in case of the limiting situations x = 0 and y —> 1 for a = 10 and 20. The 
same quantity provided by the ^(S) model is also visualized for reference. Lines of 
constant probability density [Co.o,o,o(/3)]^ Z?"* drawn as function of % and P/Pw show 
its behaviour in between the above mentioned limiting cases. 


obtained in case of which is damped by also high stiffness. Taking into acconnt 

the conclnsions drawn from fignre [H this high valne for 7 stiffness can be ascribed to 
near rotational behavionr in the gronnd band for this particular nucleus. Although the 
spectral observables of ^®°Gd are quite far from the X(3) values, it can be considered 
as quasi- 7 -rigid in virtue of the high value of the rigidity parameter. The ground and 
7 bands of all nuclei are reproduced very well, with the [3 band being the principal 
source of the discrepancy with experiment. The overestimation of the energy differences 
between consecutive fp band states is a common problem for the X(5) related models, 
which can be circumvented by considering smoothed out square well shapes for the (3 
potential [2lll25l[371|3H]. Using a polynomial potential in the (3 variable not only will 
be more natural from the geometrical symmetry perspective but will also attenuate the 
compression of the wave function against the outer potential wall which in our model is 
partially managed by increasing the rigidity y. 
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In table [2] we compare the theoretical predictions for the E2 transition rates 
computed with the values of x a from table [ 1 ] with few commonly available 
experimental data, which amount to ground-ground and 7 -ground transitions. In this 
way one can ascertain the agreement with experiment for both /S.K = 0 and /S.K = 2 
transition probabilities. In the hrst case the agreement is very good even for high 
angular momentum states, while in later case the experimental values are slightly 
underestimated. The fact that the theoretical B{E2) values vary very little between 
the considered nuclei is consistent with the presence of a stagnation interval in the /3 
transition probability density as function of y (see hgure [2]) . It is instructive to compare 
our results also with the rigid rotor estimation of the E2 transition probability 

B(B2; LK L'K’) ~ (cki-'-KK'f ■ (3.25) 

Thus, the ground-ground transitions in the present calculations are overestimated in 
respect to the rigid rotor results, while the 7 -ground ones are almost identical. The last 
result is related to the quasi-rigid description of the 7 degree of freedom. 

Another common trait of the considered nuclei is the nearly rigid rotor values of 
the experimental R 4/2 ratio which is consistent with their high experimental quadrupole 
moments [39] and microscopic mean-held calculations nnun] which also predict very 
sharp probability distributions for the /3 deformation of these nuclei. This is actually one 
of the reasons why these nuclei were found to be the best representatives of the proposed 
model. Indeed, the approximate /9-rigid behaviour of the model due to the specihc 
choice of the separated potential is in agreement with the experimental observations 
and the microscopic predictions regarding the ground state properties of these nuclei. 
Another aspect which pleads in favor of the present description are the high values 
obtained for the stiffness parameter for all three nuclei, and especially for those closer 
to rigidity, which also validate the small angles approximation made here [T8l|22]. The 
considered nuclei present however deviations from the rotational behaviour in the ground 
band when going to higher spin states. Moreover, the /? and 7 bands are not even 
remotely close to each other as it happens in the SU{3) case where these bands are 
degenerated. The experimental E2 transition probabilities for these three nuclei also 
deviate from the rigid rotor estimation as can be seen from table [2l The microscopic 
calculations of ffll also evidenced the so called centrifugal stretching phenomenon, i.e. 
the increase of the average deformation with the angular momentum. As a mater of 
fact, our model is able to describe this property due to the presence of the centrifugal 
term of the potential coming from the 7 variable. This contribution shifts the inner wall 
of the /5 square well potential, which in combination with the usual centrifugal term 
is responsible for the centrifugal stretching. A simplihed picture of this mechanism is 
described in [12] with the conhned /9-soft model [13|. It must be mentioned here that the 
inhnite square well shape adopted in our formalism is understood in the present context 
as a maximally anharmonic potential. Its usefulness is justihed also by the microscopic 
analysis [4T| regarding the 0 ^ excited states in the isotopic chains corresponding to the 
presently discussed nuclei. The results of m point to the increase in /3 anharmonicity 
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as function of the neutron number, such that its stronger for the more deformed nuclei. 
It is interesting that the aforementioned study suggest that the highest /3 anharmonicity 
corresponds to the nucleus ^^'^Gd which happens to be the best experimental realisation 
of the present model. 

The nuclei ^®^Dy and ^^^Er were previously reported as candidates for ES-X(5) 
model [22] along with ^®^Gd and ^^®Yb. Except ^®^Dy, the nomination was based only 
on the i? 4/2 ratio and the excited band heads. The inclusion of the 7 -rigidity factor 
signihcantly improved the agreement with experiment for ^®^Dy regarding the whole 
energy spectrum. Another important point regarding the comparison with ES-X(5) 
approach which is worth mentioning is the fact that present model’s hts for ^®^Gd offered 
approximately the same results, namely 7 ; ~ 0 and a = 8.5 with a = 0.415 against 
a(ES-X(5)) = 8.3 [22]. Numerical applications on the nuclei considered here were also 
performed using the exactly separable version of Bohr Hamiltonian with a Davidson 
potential [36], whose accord with experiment is poorer for ^®^Dy and ^®°Gd and better 
for ^®®Er. However, it must be mentioned that the experimental fd band considered 
in [36] for the ^®^Dy nucleus, is presently identihed as S-band [15]. The experimental 
fd band for ^®°Gd in our comparison is also different as per suggestion of m- A better 
description of these three nuclei was obtained with a Morse potential m, as well as 
with Davidson [9] and Kratzer [10] potentials combined with a deformation dependent 
mass term. However all these models posses an additional free parameter. 

4. Conclusions 

A simple exactly separable model was constructed by taking the kinetic energy of the 
Bohr Hamiltonian as a combination of prolate 7 -rigid and 7 -stable rot at ion-vibration 
kinetic operators. The relative weight of these two components is managed through a 
so called rigidity parameter %. It also dehnes the potential, such that when x = 0 the 
potential takes the separated form m(/3) -|- u{'j)/(d^, while for x —^ 1 the 7 potential 
disappears. The 7 and angular part of the eigenvalue problem associated to such a 
model is treated as in I2Z1, i.e. by applying a small angle approximation for the 7 
trigonometric functions and the corresponding potential. The later acquires a harmonic 
oscillator form centered on 7 = 0 point and characterized by a stiffness parameter a. As 
the scope of the present approach was to insert the rigidity degree of freedom into the 
description of critical point collective motion, an inhnite square well shape was chosen 
for the Id potential. This choice is consistent with the requirement that in a critical point 
of a shape phase transition the potential must be flat or having multiple degenerated 
minima [ 5 . With this choice, the (d part of the problem is brought to a Bessel type 
equation. As a result, the model’s energy spectrum and the E2 transition probabilities 
depend on two parameters, excepting the scale. Numerical applications showed that 
the energy spectrum at x 1 recovers the X(3) values regardless of the stiffness value. 
However, in the rest of the interval, the stiffness plays an important role in dehning the 
excited band heads which are rising proportionally. On the other hand, its effect on 
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Table 1. Theoretical results for ground, /3 and 7 bands energies normalized to the 
energy of the first excited state 2 + are compared with the available experimental data 
for ^®°Gd [H], [3^ and ^®®Er |3B]. The values in parentheses denote states with 

uncertain band assignment and therefore were not taken into account in the fitting 
procedure. The adimensional parameters x u are also given together with the 
corresponding deviation a defined by (13.241) . 
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13.03 

18.51 

18.41 

16.42 

16.28 

15.09 

14.71 

7 + 

4 

20.58 

20.38 

18.48 

18.25 

17.08 

16.64 

22.81 

22.62 

20.71 

20.48 

19.31 

18.82 

9+ 


25.11 

23.28 

22.96 

21.74 

21.23 

10 + 

12 + 

13+ 

14+ 

15+ 

16+ 

17+ 

18+ 

19+ 

20 + 

28.14 

27.84 

25.88 

25.68 

24.37 

23.87 


30.82 

28.98 

28.63 

27.18 

26.73 

34.31 

34.02 

32.52 

31.80 

30.14 

29.80 


37.46 

35.45 

35.20 

(32.94) 

33.09 


41.11 

39.00 

38.82 

(35.74) 

36.59 


44.99 

42.57 

42.65 


40.28 


49.08 

46.30 

46.69 


44.18 


53.38 

50.08 

50.94 


48.27 


57.90 

53.84 

55.39 


52.55 


62.62 


60.05 


57.03 


67.54 


64.90 


61.69 

X 

0.948 

0.269 

0.848 


a 

168.899 

10.309 

41.191 

a 

0.567 

0.845 

1.359 
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Table 2. Several commonly available experimental E2 transition probabilities for 
[H], [3^ and ^®®Er [30] are confronted with the model’s predictions. 

AK = 0 transition rates are normalized to the 2+ —>■ 0+ transition, while AK = 2 
transitions to the 2+ —>• 0+ transition, as in [3SJ|3^. The rigid rotor estimations are 
also presented for reference. 


Nucleus 

47 ^ 2 + 

67^4+ 

8+^67 

107 ^ 8 + 

127^107 

147^127 

2+^27 

2+^47 

2 +-)- 0 + 

27^07 

27^07 

27^07 

27->-o7 

27^07 

2+^07 

2+^07 

160 Gd 







1.87(12) 

0.189(29) 


1.45 

1.62 

1.74 

1.83 

1.90 

1.97 

1.441 

0.073 

162Dy 

1.42(6) 

1.48(9) 

1.70(9) 

1.72(11) 

1.62(20) 

1.62(20) 

1.78(16) 

0.137(12) 


1.45 

1.64 

1.77 

1.87 

1.96 

2.04 

1.444 

0.074 

166g^ 

1.44(6) 

1.71(10) 

1.72(8) 

1.80(9) 

1.71(10) 

1.84(23) 

1.86(14) 

0.151(10) 


1.45 

1.64 

1.77 

1.87 

1.95 

2.03 

1.445 

0.074 

Rigid 

1.43 

1.57 

1.65 

1.69 

1.72 

1.74 

1.429 

0.071 

rotor 










the i? 4/2 ratio has a saturation behaviour, tending to the rotational limit 3.33 at high 
values. Taking into account the above considerations, the spectral characteristics of the 
present model can be succinctly described as R 4./2 > 2.44, E{2'^) > 2.87 and 

E{2jj)/E{2'^) > 1. It is needless to say that in the lower limit of the rigidity at x = 0 
the energy spectrum recovers the corresponding ES-X(5) structure. 

Although there are two free parameters, the present solution is as poor in candidate 
nuclei as X(5) critical point symmetry and its exactly solvable analog ES-X(5). This 
is actually a common trait for parameter independent critical point solutions. Indeed, 
because the energy spectrum of nuclei change discretely with N and Z, the critical 
point of phase transition can very well be missed in between two nuclei. Moreover, 
the present model is based on the irrotational flow interpretation of the moments of 
inertia while it is well known [31] that the actual moments of inertia he in between 
these hydrodynamical realizations and the rigid ones. As a result, the search for model 
candidates revealed only three nuclei with good agreement with experiment and a high 
degree of mixing between 7 rigidity and stability. These are the N = 96 isotones ^®°Gd 
and ^®^Dy together with ^®®Er nucleus which have 98 neutrons. Considering the structure 
of the proposed model it is then not surprising to hnd its experimental realisations in 
the heavier and more deformed branches of isotopic chains undergoing shape phase 
transitions. The transitional region to which they belong is also known to exhibit level 
sequences of the collective spectrum with alternating parities which are understood 
as a consequence of the quadrupole-octupole coherent coupling [18|. Such that for a 
more realistic description of these nuclei the inclusion of at least octupole degrees of 
freedom is prerequisite. Although the hexadecapole deformation induce irremediable 
complications it must not be completely neglected given the open possibility for the 
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realization of octahedral discreet symmetry in even-even rare earth nuclei along with 
the tetrahedral one |19]. 

In what concerns numerical results, the hts for ^®^Dy designate it more closer to the 
7 -stable behavior, while for the other two the rigid excitations are stronger. Moreover, 
^®°Gd seems to be falling into the quasi- 7 -rigid category. Also our model £t for the ^®^Gd 
nucleus, conhrms it as a ES-X(5) representative [ 22 ] which in the present formalism 
amounts to y = 0. Regarding the electromagnetic properties of these nuclei, theoretical 
estimations show a good agreement with the available experimental B{E2) values for 
AK = 0 as well as AK = 2 transitions. In both cases, one observes a very weak 
dependence on the two parameters. 

Because of the above-mentioned aspects one may say that the hybrid formalism 
proposed in this paper unveils alternative features of the collective motion in the vicinity 
of the critical point of a spherical to deformed shape phase transition. The combining 
scheme described in the previous sections can be extended also to quadrupole excitations 
closer to the critical point by considering a different [3 potential, such as for example 
harmonic oscillator or Davidson [50] potentials. One expects that the predictive power 
of the model would increase based on a greater number of treatable nuclei with such 
potentials [2^136] . 
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